Inhibition of DNA ejection from bacteriophage by Mg+^ counterions 
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Experiments showed that MgS04 salt inhibits DNA ejection from bacteriophages non- 
monotonically. There is a MgS04 concentration, A'o, where the least DNA is ejected. We propose 
that this is the result of DNA overcharging by Mg"*"^ ions. As Mg"'"'^ concentration, A*', increases, 
DNA net charge changes from negative to positive. Nq is the concentration where DNA is neu- 
tral. For A'^ 7^ No, DNA is charged, prefers to be in solution to lower its electrostatic self-energy, 
hence more DNA is ejected. Our theory fits experimental data well. Mg"*"^ —mediated DNA— DNA 
attraction is found to be — 0.004fcsr per nucleotide. 

PACS numbers: 87.19.xb. 87.14.ek. 87.16.A- 



Most bacteriophages, or viruses that infect bacteria, 
are composed of a DNA genome coiUng inside a rigid, 
protective capsid. It is well-known that the persistence 
length Ip of DNA is about 50 nm, comparable to or even 
larger than the inner diameter of the viral capsid. The 
genome of a typical bacteriophage is about 10 microns or 
200 persistence lengths. Thus, the DNA molecule is con- 
siderable bent and strongly confined inside the viral cap- 
sid, resulting in a substantially pressurized capsid with 
internal pressure as high as 50 atm [J-l^. It has been 
suggested that this pressure is the main driving force for 
the ejection of the viral genome into the host cell when 
the capsid tail binds to the receptor in the cell mem- 
brane and subsequently opens the capsid. This idea is 
supported by various experiments both in vivo and in 
vitro 0, 0, The in vitro experiments addition- 

ally revealed possibilities of controlling the ejection of 
DNA from bacteriophages. One example is the addition 
of PEG (polyethyleneglycol) , a large molecule that is in- 
capable of penetrating the viral capsid. A finite PEG 
concentration in solution produces an apparent osmotic 
pressure on the capsid. This in turn leads to a reduction 
or even complete inhibition of the ejection of DNA. 

Since DNA is a strongly charged molecule in aque- 
ous solution, the screening condition of the solution also 
affects the ejection process. At a given external os- 
motic pressure, by varying the salinity of solution, one 
can also vary the amount of DNA ejected. Interest- 
ingly, it has been shown that monovalent counterions 
such as NaCl have negligible effect on the DNA ejection 
process[5|. In contrast, multivalent counterions such as 
Mg+^, CoHex+'^, Spd^"^ (spermidine) or Spm"*"^ (sper- 
mine) exert strong effect. One such results is shown in 
Fig. m where the solid circles represent experimental 
data for the percentage of ejected DNA from bacterio- 
phage A (at 3.5 atm external osmotic pressure) as a func- 
tion of MgS04 concentrationfl^ [llj. The three different 
colors correspond to three different sets of data. Evi- 
dently, the effect of multivalent counterions on the DNA 



ejection is non-monotonic. There is an optimal Mg"*"^ 
concentration where the least DNA genome is ejected 
from the phages. Similar qualitative behavior is observed 
for other multivalent counterions. 
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FIG. 1. Inhibition of DNA ejection depends on MgS04 con- 
centration for bacteriophage A at 3.5 atm external osmotic 
pressure. Solid circles represent experimental data [13, fTl| . 
The dashed line is a theoretical fit of our theory. 

In this paper, we focus on understanding the elec- 
trostatics involved in the inhibition of DNA ejection by 
Mg"*"^ counterions. We propose that the non-monotonic 
behavior observed in Fig. [T] is the result of Mg"'"^ ions 
inducing an effective attraction between DNA segments 
inside the capsid, and the so-called overcharging of DNA 
by multivalent counterions in free solution. The pro- 
posed Mg"'"^-mediated attraction between neighboring 
DNA segments inside the capsid is a central argument 
of this paper and needs to be clarified. It is well-known 
that Mg"''^ ions do not or only partially condense free 
DNA molecules in aqueous solution (l2. 13 1. However, we 
argue that due to the entropic confinement of the viral 
capsid, DNA are strongly bent and thermal fluctuations 
of DNA molecule is strongly suppressed compared to that 
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in free solution. It is due to this unique setup of the bac- 
teriophage, where DNA is pre-packaged by a motor pro- 
tein during virus assembly, that Mg"*"^ ions can induce 
attractions between DNA. It should be mentioned that 
Mg+^ counterions are shown experimentally to be able 
to condense DNA in another confined system: the DNA 
condensation in two dimension Therefore it is not 
suprising that Mg+^ ions can cause DNA-DNA attrac- 
tions inside the capsid (a zero-dimensional system). 

The overall electrostatics of Mg"^^ modulated DNA 
ejection from bacteriophages is following. Due to strong 
electrostatic interaction between DNA and Mg+^ coun- 
terions, the counterions condense on the DNA molecule. 
As a result, the net charge of DNA (77* per unit length) 
which is the sum of the "bare" DNA charges {rjbare = 
— le/1.7A) and the charges of condensed counterions be- 
comes smaller in magnitude than the "bare" charge. 
There are strong correlations between the condensed 
counterions at the DNA surface which cannot be de- 
scribed using standard Poisson-Boltzmann mean-field 
theory. Strongly correlated counterion theories, various 
experiments and simulations 15, [l6| have showed that 
when these strong correlations are taken into account, -q* 
is not only smaller than rjhare in magnitude but can even 
have opposite sign: this is known as the charge inversion 
phenomenon. Specifically, the degree of condensation, 
hence ry*, depends logarithmically on the concentration 
of multivalent counterions, N . As N increases from zero, 
77* becomes less negative, neutral and eventually positive. 
We propose that the multivalent counterion concentra- 
tion iVo, where DNA net charge is neutral, corresponds to 
the optimal inhibition due to Mg+^ induced DNA-DNA 
attraction inside the capsid. At lower or higher concen- 
trations, 77* is either negative or positive. As a charged 
molecule at these concentrations, DNA prefers to be in 
solution to lower its electrostatic self-energy (due to its 
geometry, the capacitance of DNA molecule is higher in 
free solution than in the bundle inside the capsid). Ac- 
cordingly, this leads to a higher percentage of ejected 
viral genome. The dashed line in Fig. [T] is a fit of our 
theoretical result to the experimental data for MgS04. 
The optimal Mg+^ concentration is shown to be iVo = 64 
mM. The Mg+^— mediated attraction between DNA dou- 
ble helices is found to be —0.004 /csT/base {ks is the 
Boltzmann constant and T is the temperature of the sys- 
tem). 

We begin by writing the total energy of the DNA 
molecule as the sum of the energy of the DNA segment 
ejected outside the capsid with length Lq and the energy 
of the DNA segment remaining inside the capsid with 
length Li = L — Lg, where L is the total length of the 
viral DNA genome: 



Etot(Lo) — Ein{Li) + Eout{Lo) 



(1) 



Because the ejected DNA segment is under no confine- 
ment, we neglect contributions from bending energy and 



approximate Eout by the electrostatic energy of a free 
DNA of the same length in solution. Treating the DNA 
molecule as a uniformly charged cylinder with radius a 
and linear charge density 77*, one obtains: 



EoutiLo) = -Lo{v*Vd) ln(l + rja), 



(2) 



where Z? = 78 is the dielectric constant of water and 
Ts is the Debye-Hiickel screening length of the solution. 
The negative sign signifies the fact that the system of the 
combined DNA and the condensed counterions is equiv- 
alent to a cylindrical capacitor under constant charging 
potential. The net linear charge density of DNA, 7;* , is a 
function of the counterion concentration N |l5|]: 
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-{ilj2Z) ln(iVo/iV)/ln(l + r,/a). 



(3) 



where 77c = DkBT/e is Manning critical charge density 
and Z is the counterion valence. The constant concentra- 
tion at which DNA is neutral, A^o, can be interpreted as 
the concentration of counterions next to the condensed 
counterion layer on the DNA surface. A simple derivation 
for 77* can be obtained by dividing the counterion popula- 
tion into two groups, a "bound" (condensed) counterion 
layer on the DNA and a "free" counterion population in 
solution. The distribution of the latter is assumed to 
obey Boltzmann statistics: 



N{r) = Nexp[~Zecf>{r)/kBT] 



(4) 



with (j>{r) being the electrostatic potential at radial dis- 
tance r from DNA central axis. Denoting Nq ~ N{a), 
one immediately gets 



{a) = -{kBT/Ze)ln{No/N). 



(5) 



On the other hand, the surface potential (f>{a) of a charged 
cylinder with charge density 77* in Debye-Hiickel approx- 
imation is given by [TtI : 



277* Ko{a/r,) 
D {a/rs)Ki{a/rs) 



(6) 



where Kq i are Bessel functions (this expression is twice 
the value given in Ref. 17| because we assume that the 
screening ion atmosphere does not penetrate the DNA 
cylinder). Eliminating (j){a) from Eqs. ([5]) and ([6]), one 
gets Eq. ©. 

Obtaining the concentration A^o using first princi- 
ple calculations is a complicated and non-trivial task 
[l^[l^. In general, it depends on the correlation between 
"bound" counterions at the DNA surface and its compe- 
tition with the counterion entropy. However, in practical 
situations, DNA is almost neutralized {\r]bare/il*\ ^ 1) 
by the counterions. Therefore, Nq can be very well as- 
sumed to be independent of N and 77* . Within the scope 
of this paper, we treat it as a phenomenological constant 
concentration whose value is obtained by fitting the re- 
sult of our theory to the experimental data. 
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The energy of the DNA segment inside the viral capsid 
comes from the bending energy of the DNA coil and the 
interaction between neighboring DNA double helices: 



Ein{Li,d) — Ebend{Li,d) + Eint{Li,d). 



(7) 



where d is the average DNA— DNA interaxial distance. 
To calculate Ebgnd, we employ the viral DNA packaging 




FIG. 2. A simplified model of bacteriophage genome packag- 
ing. The viral capsid is modeled as a rigid spherical cavity. 
The DNA inside coils co-axially inward. Neighboring DNA 
hehces form a hexagonal lattice with lattice constant d. A 
sketch for a cross section of the viral capsid is shown. 

model used previously [1, [H, [l^. In this model, the 
genome coils co-axially inward with the neighboring DNA 
helices forming a hexagonal lattice with lattice constant 
d (Fig. [2|). For a spherical capsid, this model gives: 



Ebend{Li,d) 

+Rln 



{-( 



V3d^ ^ V Stt 
i^+(3^/3L,dV8^)l/3 



[(i?2-(3\/3i,dV8^)2/3 



1/2}' 



(8) 



where R is the radius of the inner surface of the viral 
capsid. 

To calculate EintiLi,d), we notice that multivalent 
counterions can induce condensation of free DNA in so- 
lution 12 , 13 . DNA molecules in such a condensate are 



almost neutralized by the counterions and are arranged 
in a hexagonal lattice (similar to our viral DNA packag- 
ing arrangement) with an equilibrium interaxial distance 
do- The value of do depends on the both valency and the 
type of counterions used, but is typically about 2.8nm. 
As mentioned earlier we argue that, in the confinement 
of the viral capsid, Mg+2 also induces attraction between 
DNA segments. Neglecting finite size effect, we approxi- 
mate the interaction energy the viral DNA coiling inside 
the capsid as 



(9) 



where — e is the DNA-DNA attraction per unit length. 
Like the aforementioned parameter Nq, we treat e and dp 
as constant fitting parameters of our theory. In total, we 
have three fitting parameters (iVo, e, do) and three fitting 



constrains (the two coordinates of the minimum and the 
curvature of the curve Lo{N) in Fig. [1]). Thus our theory 
does not contain unnecessary degrees of freedom. 

Obviously, due to the strong confinement of the vi- 
ral capsid, the interaxial distance d between neighboring 
DNA double helices inside the capsid is smaller than the 
equilibrium distance dn inside the condensate. The ex- 
periments from Ref. [12| provided an empirical formula 
that relates the restoring force to the difference do — d. 
Integrating this restoring force with d, one obtains an ex- 
pression for the interaction energy between DNA helices 
for a given interaxial distance d: 



Eint{Li, d) = LiViFa (c^ + cd) exp 



'do 



-(c^ + cdo) - -(dj 



(10) 



where the empirical values of the constants Fq and c are 
0.5 pN/nni^ and 0.14 nni respectively. 

Equation ([1]) together with equations (H]), ©, ([5]) and 
pn)) provide the complete expression for the total energy 
of the DNA genome of our theory. For a given external 
osmotic pressure, Ilosm, and a given multivalent coun- 
terion concentration, N, the equilibrium value for the 
ejected DNA genome length L* is the length that mini- 
mizes the total free energy G{Lo) of the system, where 



G{Lo) = Etot{Lo) + UosmLoTTO^- 



(11) 



Here, LoTto^ is the volume of the ejected DNA segment 
in aqueous solution. The specific procedure is follow- 
ing. The energy i?i„(L — Lo,d) of the DNA segment 
inside the capsid is minimized with respect to d to ac- 
quire the optimal DNA-DNA interaxial distance for a 
given DNA ejected length, d*(Lo). Then, we substitute 
EtotiLo) = E,4L - Lo, d*{Lo)] + Eout{Lo) into Eq. ^ 
and optimize G{Lo) with respect to Lq to obtain the equi- 
librium ejected length L* {liosm ,N). By fitting this L* 
with experiment data we can obtain the values for the 
neutralizing counterion concentration, Nq, the Mg"'"^ — 
mediated DNA-DNA attraction, — e, and the equilibrium 
DNA-DNA distance do. The result of fitting our theo- 
retical ejected length L* to the experimental data of Ref. 
[lot is shown in Fig. [T] In this experiment, wild type 
bacteriophages A was used, so i? = 29 nm and L = 16.49 
fim Ilosm is held fixed at 3.5 atm and the Mg^^ 

counterion concentration is varied from 10 mM to 200 
mM. The fitted values are found to be A^o = 64 mM, 
e — 0.004 ksT per nucleotide base, and do = 2.73 nm. 

The strong influence of the multivalent counterions on 
the process of DNA ejection from bacteriophage appears 
in several aspects of our theory and is easily seen by set- 
ting d = do, thus neglecting the weak dependence of d on 
Li and using Eq. © for DNA-DNA interaction inside the 
capsid. Firstly, the attraction strength e appears in the 
expression for the enthalpy, Eq. ifTTj) . with the same sign 
as Ilosm. In other words, the attraction between DNA 



4 



strands inside capsid acts as an additional "effective" os- 
motic pressure preventing the ejection of DNA from bac- 
teriophage. This switch from repulsive DNA-DNA in- 
teraction for monovalent counterion to attractive DNA- 
DNA interaction for Mg+^ leads to an experimentally ob- 
served decrease in the percentage of DNA ejected from 
50% for monovalent counterions to 20% for Mg+^ counte- 
rions at optimal inhibition {N = Nq). Secondly, the elec- 
trostatic energy of the ejected DNA segment given by Eq. 
([2]) is logarithmically symmetrical around the neutraliz- 
ing concentration Nq- This is clearly demonstrated in 
Fig. [1] where the log-linear scale is used. This symmetry 
is very similar to the behavior of another system that ex- 
hibits charge inversion phenomenon, the non-monotonic 
swelling of macroion by multivalent counterions 21 1. 

It is very descriptive to compare our fitting values for 
e and Nq to those obtained for other multivalent coun- 
terions. Fitting done for the DNA condensation experi- 
ments by Spm+'' and Spd+'^ shows e to be 0.07 and 0.02 



fc^r/base respectively [12,|22|. For our case of Mg+^, a 



divalent counterion, and bacteriophage A experiment, e 
is found to be 0.004 ksT /haae. This is quite reasonable 
since Mg"*"^ is a much weaker counterion. Furthermore, 
Nq was found to be 3.2 mM for the tetravalent coun- 
terion, 11 mM for the trivalent counterion. Our fit of 
iVo =64 mM for divalent counterions again is in favor- 
able agreement with these independent fits. Note that, 
in the limit of high counterion valency (Z — ^ oo}, A'o is 



shown to vary exponentially with —Z^l"^ [l5l Il6l |. The 
large increase in Nq from 3.2 mM for tetravalent counte- 
rions to 64mM for divalent counterions is not surprising. 

The fitted value — e = — 0.004fcBr per base explains 
why Mg+^ ions cannot condense DNA in free solution. It 
corresponds to an attraction of — LlSfc^T' per persistence 
length. Since thermal fluctuation energy of a polymer is 
about ksT per persistence length, this attraction is too 
weak to overcome thermal fluctuations. It therefore can 
only partially condense free DNA in solution jT^. Only 
in the confinement of the viral capsid, can this attraction 
effect appears in the ejection process. It should be men- 
tioned that computer simulation of DNA condensation 
by idealized divalent counterions 23| does show a weak 
short-range attraction comparable to our e. 

The phenomenological constants — e and A^o depend 
very strongly on the strength of the correlations be- 
tween multivalent counterions on the DNA surface. The 
stronger the correlations, the greater the value e and 
the smaller the concentration Nq. In Ref. [l3|, MgS04 
salt induces strong inhibition effect. Due to this, A^o for 
MgS04 falls within the experimental measured concen- 
tration range and we use these data to fit our theory. 
MgCl2 induces weaker inhibition, thus A^o for MgCl2 is 
larger and apparently lies at higher value than the mea- 
sured range. More data at higher MgCl2 concentration is 
needed to obtain reliable fitting parameters for this case. 
In the future, we plan to complimentary our phenomeno- 



logical theory with a first principle calculation to under- 
stand the "microscopic" quantitative differences between 
MgS04 and MgCl2 salts. The authors of Ref. also 
used non-ideality as an explaination for these differences. 

In conlusion, in this letter, it is shown that divalent 
counterions such as Mg+^ can have strong effects on DNA 
condensation in a confined environment (such as inside 
bacteriophages capsid) similar to those of counterions 
with higher valency. This fact should to be incorporated 
in any electrostatic theories of bacteriophage packaging. 
The strength of short-range DNA-DNA attraction me- 
diated by MgS04 salt is first obtained by the authors. 
It is consistent with the known values for higher valence 
counterions. It provides a good starting point for future 
works with DNA-DNA condensation in the presence of 
divalent counterions. 
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